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Abstract 

Axially symmetric, stationary solutions of the Einstein-Maxwell equations with 
disconnected event horizon are studied by developing a method of explicit inte¬ 
gration of the corresponding boundary-value problem. This problem is reduced to 
non-leaner system of algebraic equations which gives relations between the masses, 
the angular momenta, the angular velocities, the charges, the distance parameters, 
the values of the electromagnetic field potential at the horizon and at the symmetry 
axis. A found solution of this system for the case of two charged non-rotating black 
holes shows that in general the total mass depends on the distance between black 
holes. Two-Killing reduction procedure of the Einstein-Maxwell equations is also 
discussed. 
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1 Introduction 


In this paper we study axially symmetric, stationary solutions of the Einstein- 
Maxwell equations having disconnected event horizon (N black holes). The so¬ 
lutions of this kind have been known since the early days of the general relativity 
Ref. [jl] and, from our point of view, are among the most interesting exact solutions 
of the Einstein equations. A great merit of these solutions is that they have a clear 
mechanical interpretation. It allows us to hope that the solutions of this class will 
help us to understand N-body problem which, we think, is one of the most fun¬ 
damental and difficult problem of the classical general relativity. We begin from 
description of the main physical properties of these solutions. 

In general the solutions of this family are parameterized by the mass, the charge, 
the angular momentum of each of the black hole and by the distances between them. 
Conical singularities at the non-extreme symmetry axis components prevent the fall 
of the black holes on each other. Interpreting these symmetry axis components 
as ’matter strings’ we get a natural definition of the interaction force between the 
black holes, namely, the tension of ’string’. It is exactly equal to the deficit angle 
of conical singularity Ref. 0](see also Ref. ||]). The interaction force of two non- 
charged black holes tends to the Newtonian limit as the distance parameter becomes 
large Refs. Jlf. 3 an d S oes to infinity as the distance parameter approaches the value 
for which two components of the horizon intersect Ref. ||. The generalization of 
these results for the case of charged black holes hasn’t been done yet. 

It is worth mentioning that for so-called extreme black holes (the horizon is con¬ 
tracted to points) the interaction force can be equal to zero Ref. ||. Unfortunately, 
this is the only known case when the balance between the gravity attraction force 
and the electromagnetic repulsion force takes place. 

The total mass of N Schwarzschild black holes in equilibrium doesn’t depend on 
the distances. However, it is a function of the distance parameters for the rotating 
black holes that is stipulated by the spin-spin interaction of the black holes. The 
charge ’densities’ of the black holes must also effect on the total mass since some 
part of the energy is needed for the charge redistribution in the final black hole. 
In the case of the static model we find that the mass of two charged black holes 
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depends on the distance parameter when their charge ’densities’ are different. To 
determine the masses of the black holes one must solve very complex constraint 
equations for the angular momenta and the angular velocities of the black holes 
Ref. [6!]. A partial solution of these constraints was found in Ref. Q. We obtain the 
equations for the physical parameters of the charged black holes in this paper. 

Despite of the presence of the conical singularity the gravity action on the sta¬ 
tionary N black holes solutions remains finite. It allows us to construct the thermo¬ 
dynamics for the system of N black holes. The thermodynamic properties of this 
system are known only for non-rotating black holes Ref. JtJ. The surface area of 
the horizon of two Schwarzschild black holes increases when the distance parameter 
decreases. 

All stationary, axially symmetric solutions with disconnected event horizon sat¬ 
isfy a certain boundary-value problem for the Einstein-Maxwell equations. The 
boundary conditions of this problem are the regularity conditions for the symmetry 
axis and for the event horizon Ref. ||]. The condition of asymptotic flatness is a 
part of the event horizon definition. The uniqueness and existence theorem for this 
problem was proven in Refs. [|], [|, 0 for non-charged black holes, and in Ref. fllj] 

for charged ones. 

The Einstein-Maxwell equations with two commuting symmetries belong to a 
wide class of completely integrable equations in the sense: they can be written as 
compatibility conditions for some auxiliary linear system of equations Refs ©IS 
p~5| , 0. Using this fact one can derive 2N-soliton solution. The 2N-soliton solutions 
do not necessarily include the family of the solutions studied in this paper. However, 
we have showed Ref. || that N black holes (non-charged) family is a subclass of the 
2N-soliton family of Belinskii-Zakharov Ref. JlTj. 

The 2N-soliton solutions of the Einstein-Maxwell equations can be constructed 
by a slight modification of the Belinskii-Zakharov method Ref. [16]. Different ap¬ 
proaches were applied in Refs. [13, 1A, 0. The 2N-soliton family of Refs. [0 0 
can not contain the black holes solutions (likely, the family of Ref. [jl8|] either). It 
was extended in Ref. [HI] with the use of the technique suggested in Ref. [0. A 
key problem with the 2N-soliton solution is that this solution has extra parameters 
which must be excluded from it to get N black holes solution. The only way to do 
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this is to satisfy the regularity conditions of the horizon and the symmetry axis by 
imposing some constraints on the parameters. In this paper we use these conditions 
at the beginning. 

In the end of this section we wish to emphasize that physical solutions for two 
black holes in equilibrium are those the interaction force of which is equal to zero. 

In the next section we state the boundary-value problem, which defines N 
charged rotating black holes in equilibrium strictly. 


2 Boundary-value problem 

In this section we formulate the boundary-value problem for the Einstein-Maxwell 
equations, which describes the all possible axially symmetric, stationary solutions 
having disconnected event horizon. Space-time manifold is said to be stationary and 
axial symmetric if there are two commuting Killing fields of which one is timelike, 
/F, and the other is spacelike, mk. The vector fields k ^ and m^ are generators of 
the isometry groups which are isomorphic R and SO(2) respectively. We denote 
these fields Kq = k and K\ = m either. By definition, the vector fields Ka satisfy 
the Killing equation and commute, 


Afi = 0 , [Ka, K b ] = 0 . 


( 1 ) 


Indexes denoted by the Latin capitals run 0 and 1. 

It is more convenient to state the conditions for the self-dual electromagnetic 
field Ft, 

Ft = F + i* F, *F liV = ^e^F a0 , 

than for the electromagnetic field F. Here, e^ va p is the alternating tensor and * 
is the Hodge operator. We assume that Ft is invariant under the action of the 
isometry groups, viz 

£j^Ft = 0, (2) 

where Ck a is the Lie derivative in the direction Ka, and the field circularity condi¬ 
tion holds, 

F^K^ = 0. (3) 
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Using F t we can write the Einstein-Maxwell equations in the form 


R\ = 8ttT^, n = = 0. (4) 

Let us introduce notations 

g AB = K A/j K% $ aA = (5) 

and a matrix-valued twist potential 

co'Xb = e a ^ (3 K AA V u K R g. (6) 

From the Maxwell equations and the condition (|2j) we get that 

$aA = V a $ A . (7) 

A simple corollary of the definitions (||) and (j|6) is 

^iiue/3^ A B K C K D = 9AC^fi9BD ~ 9AD^fi9BC- ( 8 ) 

Let us now recall the basic properties of the Killing fields 

V^V^KaA = R/3nvaK A i V^VaK^A = RapK A , 


which are consequence of Eq. (1) and the curvature tensor definition. Using them 
it is easy to check that 


V [ a W/3]AB + = 0. (9) 

Furthermore, from the Einstein-Maxwell equations we see that 

eafi^uKXR u ’ K 7 b = e^K'XF^F^Kl = -ie a ^ v K ,J A * F^F^Kl 

= i{&aB$f3A ~ $/3B$aA) = 2?'V[ q (4 , bV /3 ]<1>a). 

Hence, the identity (|9j) can be written in the following form 

V [q w ( 5]AB + 2iV [a(® B^/3]<& A ) = 0. (10) 


From Eq. (10) we conclude that there exists a matrix potential Y such that 


V oXab — U a AB + 2* < 1>_bV q ,<1>a. 
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The last equation gives 


e^Oa V Q T/i /j Kq K fj = e lw0 a^A Pj K ( cKh + 2 i$ B e^ a V a $A^- ( n ) 

It can be shown that the condition (||) is equivalent to the existence of a two- 
surface which is orthogonal to both Killing fields. From now on we chose a coordinate 
system adopted to the Killing fields, Kq = K\ = xo = t, x\ = tp. Then 

Rs 2 = labdx a dx b + gAB<lx A dx B , (12) 


where 7 a b and gAB don’t depend on t, cj). From the condition (||) we derive that F 1 ! 
can be represented in the form 

F l/3 = ^9 AB ^A[a K p]B- (13) 

Here g AB denotes ( g~ 1 )AB • The self-dual property of F 1 leads to the identity 

g AB e aa ^K'^4 = i<S> aC . (14) 


Restricting the equation O) to the two-surface that is orthogonal Ka and using 
Eq. (| 8 |) we obtain 


pecD * dY AB = QADdgBC - SAcdgsD + 2 ipecD^B * d®A- 


Here * is the Hodge operator with respect to 7 a b and 


-p 2 = det g, e 0 i = 1, e A B = ~£ba- 


Projection of the identity (14) to the same surface gives 


pg AB e B c * d$A = id® C - 


In the matrix notations, Eqs. (O) and (16) can be written as 


p * dY = —gedg + 2 ip * 


and 


(15) 


(16) 


(17) 


= ipeg 1 * (18) 

Here 4> is a column vector and <h* is its Hermitian conjugation; 4>* = (<l > o, $i)- 
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With the help of Eq. (0D and Eq. (j^) we can show that the Einstein-Maxwell 
equations take the form 


d{pG~ l * dG) = 0 , 


where 

G = 

Notice that 


1 g + 2 d>d>* d> N 

; G~ l = 

1 g - 1 -2 g- 1 * \ 

^ d>* 1/2 


—2&*g~ 1 4<f>*g~ 1 $ + 2 J 


G* = G , det G = - - 


(19) 


( 20 ) 


( 21 ) 


To prove (19), we note first that 

/ 

P G- l * dG = 


( 22 ) 


edY ied<& 

'y -2d >*edH + 2 id($*eg) -2id>*edd> 

where H is Kinnersley’s potential |22|], dH = dY + idg. So, we need only to check 
the validity of the following identities 

d$* A edd> = dd>* A edH = 0. 

From Eq. (|17|) and Eq. (|l 8 |) , it is easy to see that 


dH = ipeg 1 * dH 


Then, we have 


and 


dd>* A edd> = *d<T A e * = -dd>* A edd> = 0 


A edH = *d<S>* A e*dH = -d<T A edH = 0. 


(23) 


Here we use the standard property of the Hodge operator with respect to 7 a b and 
the equations (IS), (23). 


It is worth mentioning that the same equation (19) was derived in Ref. ||T4| 1 but 
the matrices G defined there and here are different. 

We now turn to discussion of the boundary conditions. From Eq. (0) one sees 
that d* dp = 0. Let p and z, where z defined by dz = *dp, be a coordinate system 
of the two-surface. Then 


7 a bdx a dx b = f(z , p){dp z + dz 2 ), dz = *dp, dp = — * dz. 
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One can prove that the function f(z,p ) is uniquely determined by the matrix G. 
We won’t write these equations since they aren’t need us in this paper. Virtually, 
the factor / is irrelevant to the problem of N black holes since we can compute the 
interaction force using the properties of the metric coefficient gn Ref. ||]. Let us 
remark that the Hodge operator is conformally invariant and independent of /. 

Define 

500 = —V, g\ i = X , <7oi = W. 


The event horizon and the axis of symmetry can be described as follows. The set of 
points with p = 0 and X = 0 is the symmetry axis while the set of points with p = 0 
and X > 0 is the event horizon. Denote through zk (k = 1,..., 2 N N-black holes) z- 
coordinates of the intersection points of the horizon and the symmetry axis. Then an 
interval = [z 2 i,Z 2 i-i] corresponds to the horizon of ith black hole. The symmetry 
axis components we denote by Tj for i = 1,..., N + 1 where T, = (z 2 i~i, z^i-i)- It 
is assumed that zq = —oo and Z 2 N +1 = +oo. 

We pass to a new coordinate system in a neighborhood of ith black hole, 


p 1 = (A 2 -m 2 )(l -p 2 ), rrii = 
Z2i + Z2i-1 


Z2i — Z2i-l 


z — 


= A p. 


For this coordinate system the horizon and the symmetry axis regularity conditions 
can be written in the form 

p 2 W(X,p) \ 


'i jv 


1 O' 
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($0 + fVfq) = + (A*-m?)$ 0 (A,/i), ^ = <hf + (1 — ^)<h(A, p), (25) 


(A 2 -m|)V(A,7i) 

. P 2 W(X,p) 


(l-p 2 )X(X,p), 


(24) 


where 4*0) ‘hi, X, V W are smooth functions nowhere equal to zero, and 


ncic equal uu iciu, anuL iq, 

are some constants. The second condition of ( [25[) is understood as the regu¬ 
larity condition for the electromagnetic field in the vicinity of the symmetry axis 
components, = Tf 1 in ith component (i = 1,..., N + 1) of the symmetry axis. 
The constant fh is the angular velocity of ith black hole. Note that one can derive 
the conditions for the imaginary part of T from the results of Ref. (8|] using the 
equation (18). 
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In addition, we require that 


X = p 2 (l + 0(l/r)), W = p 2 0(l/r 3 ), V = l + 0(l/r), 
$ 0 = O(l/r), $1 = 0(1), 


(26) 

(27) 


as r —» oo where r = yjz 2 + p 2 and the asymptotics (26) and (27) are differentiable. 
In particular, ( Li z = 0(l/r) and $i, p = 0(l/r). The function <3?i is given up to an 
additive constant which we fix requiring that 


lim $i(z,p) = = iq , lim $i(z,p) = T^ +1 = -iq, 

z^—oo z —>+oo 


(28) 


where q is the total charge of the system (see below). The last condition is compat¬ 
ible with the equations (|lS|). 

The charge of jth black hole is defined by an integral 


Qj = Qj - iPj = ~j ~ [ F^dS^, (29) 

47r •'Sj 

where Qj is the electric charge, Pj is the magnet monopole charge and Sj is a two- 
surface surrounding jth black hole. We assume that the two-surface belongs to a 
spacelike hypersurface which is invariant under the action of the rotational isometry 
and its normal orientation is chosen in the direction of the space infinity. Let Sj be 
a surface of revolution of a curve Cj. Then 

q j = l -jc d $ 1 = l -($f +1 -$f). (30) 

Taking this into account we see that the constant q in the normalization condi¬ 
tion ( |28| ) is the sum of the charges of all black holes q = J2 Qi- 

The mass and the angular momentum of the system are defined by integrals 

"~sjL L = lL mK “ dS ^ (31) 

where Soo is the space infinity which is a two-sphere. Using Stokes’ theorem and 
the Einstein-Maxwell equations we obtain 


and 


L = J2L h L i = ^ f C pv = m p , v - (32) 

i Hi 

M = J2 Mi, Mi = -±- [ M^dS^, AV = k KV - $ 0 Fl„. (33) 
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Here Hi is the horizon component. Since the horizon is a two-surface of revolution 
we get 


Li = -^(ynM - Yn(z 2 i-i)), Mi = j(Y 10 (z 2i ) - Yio(* 2 i-i)). (34) 


A simple corollary of the boundary conditions fl24|) and (|2q ) is the identity Yio ,2 + 
LljYii tZ \jj = 2 + Using it we find that 


Mi = rrii + 2HjLj + 4>f q t 


(35) 


Recall rrii = (z 2i - z 2 i- 1)/2. 

Under the gauge transformation <f>i —> + a the quantities Mj and L* change 
as Mi —f Mi — aQiqi and L, —> L* — respectively. We define the physical mass 

and the angular momentum of one black hole by formulas 

sA 


(J>A <J>A , <J>4 I 

ph _ M I r> „ »+1 A » rP h — r. I H-l ~ ’ 


Mf 1 = Mi + 

$ 


A 

i+1 


2 

+ <*> A 


Lf“ = L i + ft- 


(36) 


AT 


i—1 




i+1 1 

Notice that M? h and are real and exactly equal to the mass and the angular 
momentum of the isolated black hole, respectively. In addition, Mf h and L^ h satisfy 
Eq. (35) hence the product qi&f is a real quantity. 

In the end of this section we state the boundary conditions for Eqs. (©• Fn 
Eqs. (24), (25) we deduce that 


om 


pG+G- 1 = 


1 o -y 00 , 2 2<ufy 0 o, 2 ^ 

o -y 0 i >2 2^y 0 i,, 

0,2 j 


\ 0 id>o, 2 


-2 i<R 4 <f> r 


(37) 


as p = 0 and z € Tj while 


QipG'pG- 1 ^ 1 = 


\ 


yio, 2 

0 

-24>fyio,2 

Yn,z 

0 

-2$fy n , 2 

1101 

• <S> 

1 

0 



Mi = 


1 Qi 0 
0 1 0 
0 0 1 


(38) 


as p = 0 and z G I % and pG )Z G 1 = 0(1) as p —> 0 
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3 Constraint equations 


Not all of the parameters introduced in the previous section are independent. In this 
section we discuss a way of deriving the relations between the masses, the angular 
momenta, the angular velocities, the charges, the distance parameters, the values 
of the electromagnetic field potential at the horizon and at the symmetry axis. 

The system of equations (R1J) is a compatibility condition of the following pair 


of matrix linear differential equations [23] 


f-i 


Dl ^c-e- D2 „ = ^£l * 


UJ 2 + p 2 


where D\,D 2 are commuting differential operators, 


ur + p z 


(39) 


2 uj 2 

D,= 3 Z - -tr—zd L 


Pi _ a | 2 u ;/3 

9 <Juii J -^2 ~ u p~ r 9 9 '-'UJ i 

L0 Z + p z U0 Z + p- 


and u is a complex parameter that doesn’t depend on the coordinates. Studying 
the non-charged black holes in Ref. Q we have showed that well-known methods of 
investigation of completely integrable equations [24] can be applied to the problem 
of N black holes. Moreover, this problem occurs to be explicitly solvable provided 
that the relations between the parameters are found. We proceed discussing the 
principle steps of Ref. [6]. 

Let w be a root of the equation 


< jo 2 — 2 uo{k — z) — p 2 = 0, 


(40) 


where, now, k is an independent spectral parameter. From Eq. 

2uj 2 2 ujp 

d z iO = - 77 - 9 , OpUl = 


Lu 2 + p 2 ' 


uj 2 + p 2 


Hence, ip'(k) = ip(uj(k)) is a solution of the linear equation 

a fPGjG-'-upGjG- 1 

dzip (k) = -- — ^ 2 + p2 - rj; (. k ). 


it follows that 


(41) 


Invariance of Eq. (40) with respect to the transformation uj —> —p 2 /uj allows us 
to fix a branch of the function uj(k) through the inequality |cu| > p. Then, from 
Eq. (4C) we see that 


uj —> 2 (k — z), p —> 0, uj —> 2 (k — z), z —> oo. 
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After choosing a branch of u, we can define the monodromy matrix T(z,y), 
which, by definition, is a solution to Eq. ([II]) such that T(y,y) = I. Using the 

, T(z,y ) can be explicitly evaluated at p = 0. 


boundary conditions (37) and 
More precisely, for p = 0 and z, y £ Tj one obtains 

(l 


T(z,y) = 


Yoo(z)-Yoo(y) 

2 (k-y) 

n k-z _ • < j,A$o(zU'b)(y) 
u k-y k-y 


^A ^bo(z)-Vbo(y) 
* k-y 


2^(1E“+*^. 


0 


—z 


; $o(g)-$o(i/) 


2(k-y) 

while for p = 0 and z, y € Ii one gets 


K k-y 

1 + i$. 


k-y 


A &o(z)-$o(y) 

k-y 


( k-z , ®i(z)-$i(y) 

k-y ' r fc-j/ 


T(z,y) = Cl 


-l 


o \ 

U * vfc—« i 


V 


Ui(U-Ui(!t/) 

2(fc-y) 

. $i(z)-$i(y) 

2(k-y) 


1 

0 


K k—y i k—y 

x.H Yll(z)-Yu(y) 

1 k—y 

1 _ 

* fc-J/ 


The reduced monodromy matrix T(&) is defined by 


U,; 


T{k) = lim e + 1 (k,z)T(z,y)e-(k,y), 

y— 00 , 2 :—>+oo 1 


where 


/ 


e±(M) = 


\ 


V 


/ 


(42) 


1 0 0 
0 cu(fc, 2 ) t2?'<7 
0 0 1 

The matrix T(k) doesn’t depend on p. It is a consequence of the conditions 
([2?]) and d|^). Furthermore, as z € Ttv+i and y G Ti the monodromy matrix T(z, y) 
can be presented in the form 


T{z , y) = T(z, Z 2 n)T(z 2 n,Z 2 N-i) ■ ■ ■ T(z 1 ,y). 

Taking the last formula and the known expressions for T(z, y) at p = 0 into account 
one finds the limit (^) easily. We write the result in the following form 


T(k) = ATv+iTvAATv-! ... K 2 T x K\. 


( 43 ) 


Here 


Tj = 


1 - 


2 Mi 


k-Z 2 j-l 

2 tf h -i\qj\ a 

(■ k-Z2j)(k-Z2j-l ) 

Si 

k-Z2j-l 




k-z 2 j-i | 
k—Z2j ' 


k—z 2 j 

—2 qjtlj 


AMj(j>2j-l 


2(2Lf-i|<y| a )n J 2iqj 


k-z 2 j -1 

_ 2(2Lf-i| gj | a )^-i 

(k-z 2 j)(k—z 2 j-i) 


1 + 


2gj 02j - 1 
k-z 2 j-i 
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and 


/ 


Kj = 


\ 


1 Dj 0 

0 1 0 

^ 0 —i((j)2j-l — 4>2j-2) 1 ) 

where Dj = Y 00 (z 2 j-i) — Too( 3 y- 2 ) and = 4> 0 (zfc)- Recall z 0 = —oo and z 2 n+\ = 
+oo. We note also that 


<hj = ‘If - VjQf+x, Olj I = <l-f - 
A - = 2Lf 

The equations (|39| ) are invariant with respect to the transformations 


and 2Lj + <&f + iqj = 2- z^- 


Aw 




-i 


where 


/ 


= 


LO 


2|4> 0 | 2 2d>o<hi - iuj <f>o 


\ 


V 


2$i$ 0 + zw 2|«T> 1 dq 

1*0 1/2 / 
Thus, the reduced monodromy matrix T(k) has to satisfy 

/ 


T(fc) = 


V 


-/ 0 
0 i 


r*(fc) 


-/ o 
0 2 


and 


T(fe) = 


(Jy 0 


V 


o 


T*(k) 


-i 


V 


a y 0 
0 2 


(44) 


(45) 


Here I is the unit 2x2 matrix, a y is the Pauli matrix and T* is the Hermitian 
conjugation of T. The matrix T(k) can not satisfy Eq. ( fl4| ) if the physical parameters 
are arbitrary. This equation ([w|) defines the relations between the parameters. 
Unfortunately, these relations are quite complicated and we can study them in 
particular cases only. 


Let us consider one black hole first. From Eq. (44) we obtain 


Dl = 2si lMt - M24L, d 2 = asjjM, + 


u 


and 


4 >i = (4> H - 4>2 = + ^i)<7i, 


H 


iH 


L i 

Li + MiQi\qi\ 2 

2L\M\ 
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Then, the equations fl44|) and ( |35| ) are eventually reduced to the familiar formulas 
for the mass and the angular velocity of one black hole, viz 


fli = 


7i j2 2 | 2 | | 12 

Mi — vn^ + cl + |^i| , 




a = 


(Mi + mi^ + a 2 ’ ‘” i '" i ' “ ' 1,1X1 ’ “ M\ 

We can solve Eqs. (|lj) for the case of two non-rotating black holes(fli = = 0) 

with the charges of the magnet monopoles equal to zero (Pi = P 2 = 0) as well. In 
this case, L^ 11 = L ^ 1 = 0, Mf h = M\ and M 2 11 = M 2 , while Eqs. (|44|) and 
eventually reduced to a system of two equations 


are 


-Ml — 777,1 T <Pl Ql, M 2 — 7772 + <f>2 Q 21 


(46) 


where 


4>i = 


% 


1 


Ml + 7771 

1 


1 + 
1 - 


2 MiQ 2 ~ 2M2Q1 
Qi(M\ + M 2 + R) 
2M1Q2 — 2M2Q1 


M 2 + m 2 + M 2 + R ), 

and R is a distance parameter chosen such a way that zi = —mi,Z 2 = m-i,^ = 

R — 7772 and Z 4 = R + m 2 ■ For completeness, we note that 


Di = -i$f 2 , 


D 2 = 2 ), D 3 = 


and 


= 4>i Qii *% = $Q2. 


The system of Eqs. (46) is one of the main results of this paper. It gives the 
relations between parameters Mi,M 2 and 7771 , 7772 . It is preferred to consider rrij, 
which are the irreducible masses, as independent parameters of the problem. First, 
because they describe the event horizon. Secondly, in the general relativity, there is 
no notion of the energy density hence, we think, there is no sense to fix the mass of 
ith black hole. Moreover, the total mass M = Mi + M 2 contains the energy of the 
electromagnetic field surrounding the black holes as well. Let mi be independent 
parameters then in the general case the masses become functions of the distance 


parameter. These functions must be derived from Eqs. (44) and (35). 

In the particular case of Eqs. (46) we see that Mi and M 2 don’t depend on R 
only if M 1 Q 2 = M 2 Qi. This is true if 7771 Q 2 = i- The quantities Qi/ 2 nii have 
the sense of charge per unit length. Hence, if the charge density of two black holes 
treated as one dimensional objects, namely, the intervals R are different the total 
mass of the system depends on the distance between the black holes. 
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